Coherent rotations of a single spin-based qubit in a single quantum dot at fixed 
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Coherent rotations of single spin-based qubits may be accomplished electrically at fixed Zeeman 
energy with a qubit defined solely within a single electrostatically-defined quantum dot; the g-factor 
and the external magnetic field are kept constant. All that is required to be varied are the voltages 
on metallic gates which effectively change the shape of the elliptic quantum dot. The pseudospin-1/2 
qubit is constructed from the two-dimensional S — 1/2, Sz = —1/2 subspace of three interacting 
electrons in a two-dimensional potential well. Rotations are created by altering the direction of 
the pseudomagnetic field through changes in the shape of the confinement potential. By deriving 
an exact analytic solution to the long-range Coulomb interaction matrix elements, we calculate 
explicitly the range of magnitudes and directions the pseudomagnetic field can take. Numerical 
estimates are given for GaAs. 

PACS numbers: 73.21. La, 03.67.Lx, 85.35.Be 



I. INTRODUCTION 

The issue of real-time coherent control of individual 
quantum states is central to quantum computing and to 
many other ideas in the burgeoning field of quantum na- 
noelectronics. The main challenge is to isolate a system 
(or the interesting parts of a system) from its environ- 
ment in order to prevent decoherence, yet have it envi- 
ronmentally coupled enough in order to perform measure- 
ments to determine what state (or distribution of states) 
the system is in. In a solid state system — semiconductors 
in particular — encoding information in the spin, rather 
than the charge, of an electron is a promising path since 
spin couples more weakly to the environment than does 
charge. But precisely because of this weaker environmen- 
tal coupling, controlling (and measuring) the dynamics 
through external fields is slower and more problematic 
than in charge systems. 

For quantum computing, employing the spin as the 
basic qubit, for essentially the reasons mentioned above, 
was recognized early oni Here, the two-qubit gates 
are controlled electrically,^ but single qubit rotations — a 
necessary ingredient in universal quantum computing — 
require local fields (or, more precisely, local Zeeman tun- 
ing), and necessitates breaking the spin symmetry explic- 
itly. In contrast, one can define coded qubits rather 
than defining a logical qubit as being a single electron 
(or excess electron) in a single quantum dot, a single log- 
ical qubit may be defined, for example, as several quan- 
tum dots. Explicit gate sequences^ for three electrons re- 
spectively confined to three quantum dots explicitly show 
that the exchange interaction, controlled through gates 
{i.e., electrical means) alone is sufficient. This requires 
both additional gates and an order-of magnitude increase 
in gate operations. 

In the present paper, we show how a spin-based qubit, 
defined in a single quantum dot, may be manipulated ex- 
clusively by pulsing voltages applied to gates; the exter- 



nal magnetic field and the g-factor are uniform, isotropic, 
and static. Thus, both single- and double-qubit gates can 
be constructed solely through voltage pulsing with a ho- 
mogeneous, static Zeeman energy. 



II. SUMMARY 

Our qubit is encoded in the two-dimensional S = 1/2, 
Sz — —1/2 subspace of three interacting electrons in a 
two-dimensional potential well. Rotations are created by 
tuning the eccentricity of the elliptic confinement poten- 
tial. 

Any two-level system can be described as a pseudospin- 
1/2 object in a pseudomagnetic field with a Hamiltonian 
written as 



qubit 



by(Ty 



(1) 



(The most general Hamiltonian will have an additional 
term proportional to the identity operator.) The &i are 
the Pauli spin matrices, and the bi are parameters depen- 
dent upon the details of the problem. To rotate qubits, 
at least one of the three pseudofield components must be 
tunable; in principle, this degree of control can be arbi- 
trarily small. As shown below, the pseudofield for the 
present system lies in a plane, which we take to be the 
x-z plane {by = 0). In particular, we consider pseudofield 
switching between two values, bo and 6i, which differ in 
magnitude and in direction 6. 

The crucial point demonstrated below is that the 
Hamiltonian of Eq. may be realized in a single el- 
liptic quantum dot, where bx, by, and bz all have a dif- 
ferent functional dependence on the eccentricity of the 
quantum dot.^ Since this eccentricity is tunable by ex- 
ternal gates, ^ the spin-based qubit may be rotated solely 
through external gate potentials which are local to the 
quantum dot. 
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Although our results below are for two-dimensional el- 
liptic confinement, the general scheme holds equally well 
for any anisotropic (non-circular) confinement potential. 
The general requirements are guided by three considera- 
tions. First, the two qubit states |0) and |1) should both 
have the same spin (1/2) and spin projection. Second, if 
the two states differ by at least one spin- flipped pair, the 
relaxation should then be governed by the spin (rather 
than charge) relaxation time, regardless of the orbital 
configurations. Third, if those spin-1/2 states which de- 
fine the qubit are the two lowest-energy states, then one 
can serve as the initial state, prepared by equilibration. 

In the following section, we outline an exact solution to 
the one-body problem. This solution has been published 
before,® but we provide an alternate derivation based on 
Bose operators, similar to the circular case, which will 
facilitate the second-quantized treatment with interac- 
tions. 

In section IIVI we consider interactions. We provide, 
for the first time, an exact, closed- form expression for all 
Coulomb matrix elements (in the single-particle eigenba- 
sis), valid for arbitrary quantum numbers. 

We next detail the explicit construction of our qubit in 
section Ivl and derive Eq. giving expressions for the 
pseudofields in terms of the various exchange energies, 
and, ultimately, in terms of the parameters appearing in 
the electronic Hamiltonian. 

Following this, we give an explicit sequence of confine- 
ment deformations which enables a qubit flip and give 
estimates based on GaAs lateral dots using realistic po- 
tential and material parameters. 

III. ONE-BODY HAMILTONIAN: EXACT 
SOLUTION 



from which the adjoint operators (aj, 0,2) can easily be 
found. These four operators satisfy the canonical Boson 
commutation relations. The dimensionless parameters 
a±, (3± are defined by 

/?.-(l±^j ' (4b) 

and we have also defined the (hybrid) magnetic length 
= hl{mujf)), cyclotron frequency uJc = eB/{mc), as 
well as^ 

uo = K + 2 K + co'y)] "'^ , = H - u^l) , (5) 

n = {ut + u:lu:lf\ (6) 

The Bose operators of Eq. © diagonalize the elliptic 
Hamiltonian, Eq. 

H = nn+{a\ai + \] +nn^[ala2 + ]^ , (7) 



where f2± = ^ y/uj^+u}^±2if^ . (In the isotropic limit 
of ujx LUy, we have a± — > 1 and f3± — > 1. The Bose 
operators and the Hamiltonian then reduce to the usual 
isotropic onesiiS) 



The Hamiltonian for a noninteracting elliptic quantum 
dot is given by 

We have neglected the Zeenian term since it plays no 
significant role in what follows. Equation (0) describes 
one electron trapped in a plane, under a perpendicu- 
lar magnetic field — we use the symmetric gauge, A = 
B{—y,x,0)/2 — with further lateral confinement by two 
different parabolic potentials with frequencies lu^ and LUy. 
This describes an elliptic confinement with the rotational 
symmetry (and consequent angular-momentum conserva- 
tion) explicitly broken. 

Equation jSJ may be diagonalized by introducing Bose 
operators analogous to the isotropic case. (For an al- 
ternative but equivalent solution to the elliptic one-body 
problem, see Ref.i^). These operators are explicitly given 
by 
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l/a- 



x/2£o 
Pyio/^ 



Y 



(3a) 
(3b) 



IV. COULOMB MATRIX ELEMENTS: 
SOLUTION 



EXACT 



For the electron interactions, we use the long-range 
Coulomb energy (~ 1/r) and work in the second quan- 
tized formalism using the exact single-particle basis \mn) 
{a\ai\mn) = n\mn), d\d2\mn) = m\mn)); hence Vc = 
5 1] Vjjfcic|^c]^,c/cr'Cfecr, where all indices (ijklacr') are 
summed over; each Latin index represents a pair of or- 
bital quantum numbers (to, n) and the Greek indices rep- 
resent spin (cr, a' — ±1/2). Calculation of the matrix ele- 
ment Vijki proceeds through the two-dimensional Fourier 
transformpi 



V^jki ^ y d^g^— (mini,TO2n2|e*'^'(''i ''^^ Itos^s, TO4n4), 

(8) 

by writing the position operator r — {x,y) in terms of 
the Bose operators in Eq. Q and their adjoint. After 
some calculation, we obtain 
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2 Ur, n \ min(ni,n3) / \ / \ min(r?ii,m3) , , , 

e iKl-niQ) sr^ I ( "1 W "-3 \ , / mi \ ( 



2 11/0=1 ™fe''^fc' ^1=" \ / \ / p2^o ^ ^ ^ 



min(n2,n4) / \ / \ min(m2,m4) / \ / \ /■! \ i \ 



, (9) 



P3=0 \ / \ / P4=0 



where 9^ = 5 + — 2pi). The integral 

may be expressed as a sum of elementary functions and 
complete elliptic integrals of the first, second, and third 
kinds. The function A = \(x) is explicitly given by 



(|u|2 + |^;|2)9-' + ^ 



(10) 



where n^j = ^^ + ^^-(^1+^3), ■m^j = rm + mj - {P2+P4), 
and 



a_/3_ 



(11) 



The matrix element, Eq. (|SJ), vanishes if J^ii^i + '^0 
is odd and is real otherwise. In the isotropic limiti 
(lUx — LOy) the expression simplifies considerably and 
conservation of angular momentum emerges explicitly. 
Equation ^ is an exact result, valid for any set of quan- 
tum numbers m^, (i = 1, . . . , 4). It can be used as the 
basis of a numerical treatment of the many-body prob- 
lem. 



V. QUBIT CONSTRUCTION 

Rotations are enabled through the mutual exchange 
interactions among the confined electrons. In what fol- 
lows, we consider three-particle antisymmetric state vec- 
tors of the form |mini(Ti, m2n2(T2, manscrs), with fixed 
orbital states (m^, rii). For a given set of orbital quantum 
numbers, we construct the qubits from the (exact) two- 
dimensional subspace of the three-electron problem with 
spin 5 — 1/2, Sz = —1/2. We shall consider the three 
orbital states (m, n) = (0, 0), (1, 0), (2, 0) with no double 
occupancy. We stress, however, that neither single occu- 
pancy nor three orbital states (only) are essential to the 
main conclusions. The important point is that the spin- 
degenerate space is two-dimensional — an exact result — 
and that the shape of the dot is tunable — an experimen- 
tally demonstrated fact.^ The resulting eight-dimensional 
Hilbert space is spanned by the antisymmetrised (Slater 
determinant) states |00(To, IOcti, 20cr2), which we will sim- 
ply write as |tTo: ""i, 172) (but note that these are antisym- 
metrised states). Three spin-1/2 particles can be com- 
bined to form a spin-3/2 quartet and two orthogonal spin- 
1/2 doublets. The two \S,Sz) = |l/2,-l/2) states are 



orthogonal and form our two qubit states |0) and |1) 
They are explicitly given by 



|0)^— (2|UT)-Im)-|TU»: 



ID 



V2 



(UTi>-|Tii». 



(12a) 
(12b) 



These states are linear combinations of single- 
determinant state vectors and, as such, go beyond 
the standard Hartree-Fock treatment. What's more, 
at finite magnetic field, these states are both lower in 
energy than spin-1/2 states involving a doubly occupied 
s-shell. 

We project the total Hamiltonian — consisting of both 
one-body, Eq. {T)), and two-body, Eq. ||HJl, terms — down 
to our two-dimensional qubit subspace, spanned by the 
vectors |0) and |1). This can be mapped to a pseudospin- 
1/2 problem whose general form is given by Eq. The 
pseudomagnetic field components are given by various 
exchange interactions. We find hy = 0, whereas 

= — (V022O - V1221) , (13a) 

bz = -Voiio + ^ (V1221 + Vb22o) • (13b) 

The pseudofields hx and bz depend on different combina- 
tions of exchange- interaction matrix elements, and each 
of these depends differently on the ratio r = ujy/ujx- This 
will be true of almost any anisotropic confinement poten- 
tial. Because of this, the direction of the pseudofield can 
be changed — inducing coherent rotations of the qubit — 
by changing the anisotropy parameter r. Analytic ex- 
pression for the various exchange energies in Eq. (|13|l are 
given in the Appendix. 

Figure^shows the angle 9 of the pseudofield h (relative 
to the positive x axis) as a function of both anisotropy r 
and (actual) magnetic field z = uJc/oJx- The larger values 
of r are the physically relevant ones. (The isotropic case, 
LOx = ^y corresponds to r = 1, whereas r = is the 
one-dimensional limit.) The figure shows that at a fixed 
magnetic field z, a range of pseudofield directions are 
available for qubit rotations by varying the voltage-tuned 
anisotropy r. In both extremes, r = 0,1, Fig. ^ shows 
no dependence on 9 with magnetic field z\ in both cases, 
the system essentially has only one tunable parameter 
which, in the logical qubit space, tunes the magnitude 
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Pseudofield angle G [degrees] Pseudofield magnitude b [e^(2/7T)^^^/(32lo)] 




O 0.2 0.4 0.6 0.8 1 O 0.2 0.4 0.6 0.8 

Z — COj, / COj^ z — CO^ / (Oj^ 



FIG. 1: Contour plot showing the angle 6 (left plot) and the magnitude b (right plot) of the pseudomagnetic field 6 as a function 
of quantum dot anisotropy r — Uy/uix and (actual) magnetic field z — ujc/lox- The angle 9 is measured from the positive x 
axis and b lies in the x-z plane. Angles are measured in degrees and magnitudes in units of e^y'2/7r/(32i?o). Shown are lines 
of constant 6 and b respectively. 



of the pseudofield (through the hybrid magnetic length 
£q). Figure ^ also shows how the magnitude (in units of 
^ 2 1 -K I {'yilo)) of the pseudofield changes as a function 
of r and z. In general, both the magnitude and direction 
of the pseudofield are altered by the anisotropy. 



VI. EXPLICIT QUBIT FLIP SEQUENCE 

By tuning r{€) in real time, a qubit flip |0) — > |1) can be 
performed; we give here an explicit example. It is useful 
to rotate our qubit, Eq. 112|l . so that it is oriented parallel 
(and antiparallel) to the direction of the pseudofield for 
r = 1, given explicitly by 

Thus, our rotated qubit states are |0) = c_|0) — c+|l) 
and |i) = c+|0) + c_|l), where c± = ^(6o ± 6oz)/(26o). 

The initial (t = 0) qubit state is along the pseudofield 
direction bg given by ro — 1, which, in our rotated frame, 
we take to lie along the z axis. The field is then pulsedi^ 
to a new value hi given by, ri < 1. (This field lies in the 
x-z plane.) The qubit will precess about bi with period 
Ti = Trh/bi. Half a period later, a,t t = Ti/2, the qubit is 
again in the x-z plane, whereupon the field is pulsed back 
to 6o- The qubit precess about this new field with period 
To = TTh/bo- Half a period later, at t = (Ti + To)/2, the 
field is again pulsed to bi and the process is repeated 
every half period. (Actually, the pseudofield does not 
need to be switched every half period; an odd number of 
half-periods suffices.) If the angle fi between bo and bi is 



chosen such that /i = t: /{2k), where k is an integer, the 
qubit may be flipped by k pulses at bi with pulse width 
Ti/2, each separated by an interval Tq/2 at bp. The total 
switching time is t^"^ = kTi/2 -\- {k — \)Tq/2 and can be 
very fast. (See below). The qubit can in fact cover the 
Bloch sphere by judicious choice of pseudofields, which 
are entirely controlled by the quantum dot anisotropy. 

For definiteness, we give here numerical estimates 
based on material parameters for GaAs. We take lOx = 
6 meV, while uiy switches between 3 and 6 meV. We also 
take a (fixed, uniform) magnetic field of _B = 0.42 T. 
Thus, r = 1, 0.5 and z = 0.12 for GaAs. For r = 1, 
the pseudofield is explicitly given by Eq. I|14|l and yields 
a magnitude of 5o ~ 1.61 meV. At r = 0.5 the magni- 
tude is decreased, hi w 0.94 meV, whereas the direction 
is increased. Figure |21 shows both the direction and 
magnitude of the pseudofield for these particular param- 
eters. The field bi is tilted away from bo by = 9°. 
This gives a qubit flip in ten pulses. With these pseudo- 
field values, the precession periods are Tq = 1.3 ps for bo 
and Ti = 2.2 ps for bi. The lower hound on the flipping 
time i^'P is for a pseudofleld switch every half-period; this 
yields i^'P ~ 16.8 ps. These times are closer to optical 
frequencies than what is currently achievable using pulse 
generators. Recent pulsed-gate experimentsi^ employed 
electrical pulse-widths on the order of 10 ns. With such 
pulse generators, we have t^'P « 190 ns. 



VII. DISCUSSION 

Our qubit, Eq. (|12|l . has been constructed from a linear 
combination of single-determinant (Hartree-Fock) state 
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FIG. 2: Direction and magnitude of magnetic field {by = 0) 
for various anisotropy values r — uiy /ujx sd, fixed magnetic field 
LJc = 0.12tJa;. This is essentially a plot of Eq. 11311 in units of 
\J 2 1 -K I (ZUq) . The hatched circle and ellipse are schematics 
of the quantum-dot shape at different anisotropics. 



vectors, where the orbital degrees of freedom have been 
frozen out. But the general scheme is certainly not lim- 
ited to our specific state vectors. In general, each logical 
qubit state can be written as a correlated many-body 

state |Q> = Eitt^'^V)*- where \Q) ^ |0), |1) is the logi- 
cal qubit state and the are antisymmetrised orthonor- 
mal states, j?/;) = jmiriicri, TO2n2C2, niaTiscra), such that 
\Q) is a spin eigenstate with S = 1/2, Sz ~ ^1/2. 
Equation H12|) . for example, has (mini, m2n2, Wsna) = 
(00,10,20) for both \Q) = |0), and \Q) = |1); the dif- 
ferences between the two logical states are, in this case, 
solely due to spin flips and phase factors of ±1. Although 
there is no requirement that the orbital degrees of free- 
dom are identical for each qubit state, it is neverthe- 
less advantageous to have the orbital quantum numbers 
identical since this will reduce the electromagnetic fluc- 
tuations which would be present if the qubit rotation 
involved orbital transitions as well as spin transitions. 

It is always possible to deflne the logical qubit states 
in such a way that they differ only by spin flips and rel- 
ative phases and not by their orbital quantum numbers. 
This statement is not restricted to the simple (yet rele- 
vant) case of that described by Eq. (|12|l . It is an exact 
result, valid even for correlated states involving many 
Slater determinants. Thus, voltage fluctuations due to 
orbital transitions can be mitigated. 

It is also possible to choose the qubit states such that 
one is the ground spin-1/2 state and, consequently, state 
preparation can be a matter of equilibration. 

Finally, the two qubit states will not be energetically 
degenerate. Thus, each qubit state will have differ- 
ent transport characteristics; the magnitude of current 
through the dot will depend differently on gate and bias 
voltages for each of the qubit states. This may be ex- 
ploited to be used as a detection scheme for flnal readout. 
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APPENDIX: EXCHANGE ENERGIES 

The pseudofleld, Eq. H13|l . is determined by various 
exchange energies. These Vijki are in turn determined 
from the exact expression of Eq. © with the subscripts 
{i,j,k,l) = (mi, TO2, ma, TO4) and all = 0. For the 
cases on interest here, the relevant Vijki are given by: 

Foiio = CX2, ^0220 = \cXi, (A.la) 

1^1221 = - 2X4 + 4X2^ , (A.lb) 

where C = e^/(47r^o) is the Coulomb energy scale, Xg = 
2(^+3)/2r((s + l)/2)/,, and 



h = du 



(l-w2)l/2(au2^_^)(.-Hl)/2- 



(A.2) 



Each Is is a linear combination of complete elliptic inte- 
grals of the flrst and second kindi^ 



Is = AsK{m) + BsE{m), 



(A.3) 



where m = {at ~ at ) / , and the coefficients As and 



+ 

Bs are given by 



A, = 



1 



(A.4a) 
(A.4b) 



a^VbV 3b{a + b)^ 
_ I fs ciy^ Aiy^a + 2b) \ 
'^a^VbV b{a + b) 1562(a + 6)2j' ^""'^'^ 



Bo = 



aVb{a + b) 



( v{a + 26) 
= —i^, I c ' 



a2V6(a + &) V 36(a + fe) 
1 (ic^v 2c//2(a + 2&) 



a^-JbKa + b b{a + by 

8v^{a + 2by 3i/3 



15b^{a + by 5b{a + by 



(A.4d) 
(A.4e) 

, (A.4f) 



where 



a — 


w+ 




6 = 


Pl/al - 


^alpt, 


c = 


atP'_ - 


0Xla\, 


d = 






V — 


ad — be, 





4) 



(A.5b) 
(A.5c) 
(A.5d) 
(A.5e) 
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and the a± and P± are given in Eq. 



URL: http : // soliton . phys . dal . ca 

D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998), cond-mat/9701055. 

G. Burkard, D. Loss, and D. P. DiVincenzo, Phys. Rev. B 
59, 2070 (1999), cond-mat/9808026. 

D. Bacon, J. Kempe, D. A. Lidar, and K. B. Whaley, Phys. 
Rev. Lett. 85, 1758 (2000). 

J. Kempe, D. Bacon, D. A. Lidar, and K. B. Whaley, Phys. 
Rev. A 63, 042307 (2001). 

D. P. Divincenzo, D. Bacon, J. Kempe, G. Burkard, and 

K. B. Whaley, Nature 408, 339 (2000). 

In the istropic limit of a circular quantum dot, we find 

by = and bx (x bz; only the magnitude, not the direction, 

of the pseudofield can be changed; the direction of the 

pseudofield can only be changed through control of the 

anisotropy. 

J. Kyriakidis, M. Pioro-Ladriere, M. Ciorga, A. S. Sachra- 
jda, and P. Hawrylak, Phys. Rev. B 66, 035320 (2002). 
A. V. Madhav and T. Chakraborty, Phys. Rev. B 49, 8163 
(1994). 

Without loss of generality, we restrict lux > i-Oy in what 



follows. 

L. Jacak, P. Hawrylak, and A. Wojs, Quantum Dots 
(Springer, Berlin, 1997). 

The notation | • • • ) denotes unsymmetrised states, while 
I • • • ) denotes properly antisymmetrised and normalised 
states. 

A smoother pulse shape for r{t), while more realistic, 
would entail a careful treatment of time-ordering and 
would unnecessarily obfuscate the present simple example. 
T. Fujisawa, D. G. Austing, Y. Tokura, Y. Hirayama, and 
S. Tarucha, J. Phys.: Condens. Matter 15, R1395 (2003). 
We define elliptic integrals of the first and second kind as 



K(m) = ( dt ^ 
Jo \/(l-t^)l 



V(l-t^)(l-mt2) 



and 



E{m) = dt 



I 



^ ,^ VI - mt2 



respectively. 



